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ABSTRACT 

We study the observational constraints on the growth of massive black holes (BHs) in 
galactic nuclei. We use the velocity dispersions of early-type galaxies obtained by the 
Sloan Digital Sky Survey and the relation between BH mass and velocity dispersion to 
estimate the local BH mass density to be p.(z = 0) ~ (2.5±0.4) x 10 5 /iq 65 M© Mpc -3 . 
We also use the QSO luminosity function from the 2dF Redshift Survey to estimate 
the BH mass density accreted during optically bright QSO phases. The local BH mass 
density is consistent with the density accreted during optically bright QSO phases if 
QSOs have a mass-to-energy conversion efficiency e ~ 0.1. By studying the continuity 
equation for the BH mass distribution, including the effect of BH mergers, we find 
relations between the local BH mass function and the QSO luminosity function. If the 
BH mass is assumed to be conserved during BH mergers, comparison of the predicted 
relations with the observations suggests that luminous QSOs (iboi > 10 46 ergs _1 ) 
have a high efficiency (e.g. e ~0.2, which is possible for thin-disk accretion onto a 
Kerr BH) and the growth of high-mass BHs (> 10 8 M ) comes mainly from accretion 
during optically bright QSO phases, or that luminous QSOs have a super-Eddington 
luminosity. If luminous QSOs are not accreting with super-Eddington luminosities and 
the growth of low-mass BHs also occurs mainly during optically bright QSO phases, 
less luminous QSOs must accrete with a low efficiency < 0.1; alternatively, they may 
accrete with high efficiency, but a significant fraction should be obscured. We estimate 
that the mean lifetime of luminous QSOs (L bol > lO^ergs" 1 ) is (3-13)xl0 7 yr, which 
is comparable to the Salpeter time. We also investigate the case in which total BH 
mass decreases during BH mergers due to gravitational radiation; in the extreme case 
in which total BH entropy is conserved, the observations again suggest that BHs in 
most luminous QSOs are Kerr BHs accreting with an efficiency > 0.1. 

Key words: black hole physics - galaxies: active - galaxies: evolution - galaxies: 
nuclei - quasars: general - cosmology: miscellaneous 



1 INTRODUCTION 



Most nearby galaxies contain massive dark objects at their centers (e.g. Kormendy & Richstone 1995, Kormendy & Gebhardt 
2001), which are presumably black holes (BHs). The existence of these objects was predicted by arguments based on quasi- 
stellar object (QSO) energetics and demography (e.g. Soltan 1982, Rees 1984). Studies of central BHs in nearby galaxies have 
also revealed a tight correlation between BH mass and galactic velocity dispersion (Fcrrarese & Merritt 2000, |Gcbhardt et 



al. 200C ) , and a less tight correlation between BH mass and the luminosity (or mass) of the hot stellar component of the host 
galaxy (e.g. Kormendy & Gebhardt 2001 and references therein; by "hot" component we mean either an elliptical galaxy or 
the bulge of a spiral or SO galaxy). These correlations strongly suggest a close link between the formation and evolution of 

galaxies and their central BHs. 

The simplest and most elegant argument relating the properties of QSOs to BHs in nearby galaxies is due to [goltan 



(1982), who pointed out that the luminosity function of QSOs as a function of redshift traces the accretion history of these 
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BHs. The BH mass density due to bright QSO phases can be estimated directly from the local energy density in QSO photons 
and an assumed mass-to-energy conversion efficiency (Soltan 1982; Chokshi & Turner 1992). Comparison of this density with 
the local BH mass density provides considerable insight into the formation and growth of massive BHs (e.g. Haehnelt & 
Kauff nann 2001 and references therein). 

In § y, we provide new estimates of the local BH mass density and the BH mass density accreted during bright QSO 
phases. There are several reasons why a re-analysis of this classic problem is timely, (i) The past several years have seen 
dramatic improvements in optical QSO surveys, both in terms of the total numbers of QSOs and the parameter ranges (i.e., 
luminosities and redshifts) over which the QSO luminosity function is reliably determined. For example, the 2dF QSO redshift 
survey flBoyle et al. 2000[ ) and the Large Bright QSO Survey flHewett, Foltz fc Chaffee 1995[ ) have found ~ 6000 QSOs 
with absolute magnitude —26 < Mb < —23 (fl m = 1,£Ia = 0, h = 0.5), at redshifts 0.35 < z < 2.3; the Sloan Digital Sky 
Survey (SDSS) has found several hundred QSOs with redshifts 3.6 < z < 5.0 (e.g. Fan et al. 2001; see also other surveys at 



high redshifts, e.g. Schmidt, Schneider & Gunn 1995 etc.) (ii) The SDSS has provided a sample of ~ 9000 nearby early-type 
galaxies (elliptical galaxies and SO galaxies; Bernardi et al. 2001) with accurate luminosities and velocity dispersions. Using 
the tight correlation between BH mass and galactic velocity dispersion, we can estimate the BH masses in these galaxies 
and thereby obtain the mass function of BHs at z — 0. (iii) A large number of obscured active galactic nuclei (AGNs) (e.g. 
~4-10 times the number of unobscured AGNs) is required to explain the observed X-ray background spectrum (e.g. Fabian 
& Iwc sawa 1999) ; Gilli, Salvati & Hasinger 2001). The BH mass density inferred from the X-ray background by Fabian & 
Iwasava (1999) is larger than the one derived from optically bright QSOs by Chokshi & Turner (1992) by a factor of 3-4. 
BHs in local galaxies are the remnants of both obscured and unobscured AGNs. Thus, an accurate estimate of the local BH 
mass density and the accreted mass density from optically bright QSOs can also constrain the fraction of obscured AGNs for 
comparison to models of the X-ray background. 

If BH growth occurs mainly during optically bright QSO phases, then not only should the BH mass density in local 
galaxies be consistent with the total energy density in QSO light, but also the BH mass distribution in local galaxies should 



be consistent with the QSO luminosity function (e.g. small fc Blandford 1992; Ciotti, Haiman fc Ostriker 2001 



Marconi 



& Sal jrati 200 H ). In § §, by studying the continuity equation for the BH mass distribution, we test this hypothesis and give 
more general observational constraints (e.g. on the mass-to-energy conversion efficiency and lifetime of QSOs) on the growth 
of massive BHs in galactic nuclei. The effects of BH mergers (caused by galaxy mergers) on the BH mass distribution are 
also considered. Since the physics of BH mergers is poorly understood, we consider two extreme cases: one is the classical 
case (mergers do not emit gravitational radiation and the total BH mass is conserved during BH mergers), which is implicitly 
assumed in Soltan's argument; the other is the adiabatic case (total BH area and entrop y is conserved during BH mergers 
and mass is radiated away as gravitational radiation; see also Siotti fc van Albada 200l| ). Our conclusions are summarized 

In this paper, the Hubble constant is written as Ho = lOO/i km s 1 Mpc ; and if not otherwise specified, the cosmological 
model used is (f2 M] fi A , h) = (0.3, 0.7, 0.65) (|Wang et al. 200C|). 



2 TOTAL BH MASS DENSITY 
2.1 Early-type galaxies in the SDSS 

The SDSS will image ~ 10 4 square degrees of the sky in five bands (u, g,r,i, z) with central wavelengths (3560, 4680, 6180, 



7500, 8870) A, and take spectra of ~ 10 6 galaxies and ~ 10 s QSOs (York et al. 2000). Among the spectra to be taken by 
SDSS, there will be roughly 2 x 10 5 spectra of early-type galaxies. A sample of nearly 9000 of these with well-measured 
velocity dispersions (<51ogu ^0.02-0.06 dex with a median value ~ 0.03 dex), in the redshift range 0.01 ^ z ^ 0.3, has been 



selected from early SDSS observations by Bernardi et al. (2001 ). The joint distribution of luminosities, sizes and velocity 



dispersions for this sample is well described by a tri-variate Gaussian distribution in the variables M = —2.5 log L, R = log R , 
and V — logo - , where L is the total luminosity, R is the effective radius and a is the line-of-sight velocity dispersion within 
a circular aperture of radius R /8. A maximum-likelihood fit to the sample gives the mean values of the above variables, 
(M) = —19.65 + 5 log h — 1.15z (the 1.15z term represents evolution), (R) = 0.36, and (V) = 2.20 in the g* band (here 
10 H and 10 v are in units of h 1 kpc and kms , respectively; g* is used rather than g, similarly for other bands, because 
the photometric calibration is preliminary; for details, see Stoughton et al. 2002) and their standard deviations <tm = 0.84, 
/jr = 0.25 and av = 0.11 (cf. table 2 in Bernardi et al. 2001). By extrapolating this luminosity distribution at both faint and 
bright ends, the average comoving number density of early-type galaxies is found to be 

(f>* = (5.8 ± 0.3) x 10~ 3 /i 3 Mpc" 3 . (1) 

For comparison to other studies, we must convert magnitudes in the SDSS g* band to the Johnson-Morgan B band. We take 
the color transformation to be: 
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B = g* + 0.47( 5 * - r*) 



0.17, 



(2) 



which is derived from the preliminary SDSS calibration ( Smith ct al. 2001 ) of the Landolt stars ( Landolt 1992a| Jb|). The 



standard deviation in this relation is 0.03 mag. The mean color of the early-type galaxies in Bernardi et al.'s SDSS sample 
is g* — r* ~ 0.74 — 0.30z with standard deviation 0.06 mag in the colors of individual galaxies (cf. table 6 in Bernardi et al. 
2001). 



2.2 Local BH mass density 



The SDSS data provide us with the distributions of velocity dispersion and luminosity in nearby early-type galaxies. We denote 
the luminosity function and the velocity-dispersion distribution of early-type galaxies as n^ arly (L, z) and n%^ Tly (a, z), where 
n^ arly (L, z)dL [or n^ arly (cr, z)da] represents the comoving number density of early-type galaxies in the range L —* L ± AL (or 
<r — > (T + der) at redshift z. The distributions n^ lly (L, z) and n^ lly (a, z) are obtained by integrating the tri-variate distribution 
in M = — 2. 5 log L, R — \ogR and V = log a over two of these variables. Thus, the average number density of these galaxies 
is (eq.|) 



early 



L,z = 0)dL = 



early 



(a, z = 0)dcr, 



and the luminosity density is 



J = 



LnT ly {L,z = Q)dL. 



(3) 



(4) 



= 6.1 X 10 h Lg*,© Mpc . For comparison, the total 
Mpc -3 (Blanton et al. 2001). This total luminosity 



The luminosity density derived from the Bernardi et al. sample is j g 
luminosity density derived from SDSS is j 9 * = (2.8 ± 0.4) x 10 s h L g * i0 
density is a factor of 1.5 to 2 times higher than earlier estimates, but Blanton et al. argue that the difference arises because 
earlier surveys neglect low surface-brightness regions that are included in SDSS. The luminosity density in the Blanton et al. 
sample can be compared directly to the Bernardi et al. sample, as follows. Strateva et al. (2001) show that early- and late-type 
galaxies in SDSS can be separated by a cut in concentration index at C — 2.6; applying this cut to the Blanton et al. sample 
we derive j g * = 7.4 x 10 7 /i Lg*,© Mpc -3 for early-type galaxies. This result differs by only 20% from the estimate of Bernardi 
et al., which uses several parameters including the concentration index to separate early- from late- type galaxies. Thus the 
normalization of the sample isolated by Bernardi et al. appears to be co nsistent with other ear ly-type galaxy samples in SDSS. 
The relation between BH mass and galactic velocity dispersion is ( Tremaine et al. 2002 ) : 



M m = (1.5 ±0.2) x 10 8 M ( 



( °° v 

V200 kms" 1 / 



(5) 



where a± — 4.02 ± 0.32, M, is the BH mass, and a e is the luminosity-weighted line-of-sight velocity dispersion within a slit 
extending to the effective radius. Note that <t c in equation (JE)) is the velocity dispersion within a slit extending to the effective 
radius R ( jGebhardt et al. 2000| ) while a in the SDSS is the velocity dispersion within a circular aperture extending to R /8; 
however, replacing a c with a in equation (|5|) will not caus e much difference as the two definitions should give very similar 
results (Tremaine et al. 2002). Merritt & Ferrarese (2001a) give an alternative version of the correlation (||) that is based on 
the velocity dispersion within R /8, and find a steeper slope, «i = 4.72. The reasons for the difference in slopes are discussed 
by Tremaine et al. (2002). We have performed all of the calculations below using both versions of the correlation and the 
differences are negligible. 



The mean relation between the _B-band luminosity of the hot stellar component and the BH mass is (Kormendy 



Gebhardt 2001) 



M. = 0.7 x 10 s Mr. 



( L i 



\W 10 L Be 



(6) 



where = 1.08 (Kormendy & Gebhardt 2001). The scatter in relation (j6J) is substantially larger than the scatter in relation 
(^). The statistical uncertainties in 02 and in the normalization are about ±0.15-0.2 and 30%, respectively. A similar relation 
to equation (^) was derived by Magorrian et al. (1998), but with normalization larger by a factor of 3-4, presumably because 
of the limitations of ground-based spectroscopy and isotropic dynamical models. 



Most of the galaxies in Tremaine et al. (2002) and Kormendy & Gebhardt (2001) have distances derived from Tonry et 



al. (2001), which implies a Hubble constant h = 0.73±0.04±0.11 (Blakeslee et al. 2002), and the normalizations in equations 
(^) and (|^) have been adjusted to our assumed Hubble constant h = 0.65. 

Because the mass-dispersion relation is the tighter of the two relations, we will first use equati on (|e|) to estimate the local 
BH mass density. Thus, with the distribution of velocity dispersion n^ Tly (a, z — 0) determined by Bernardi et al. (2001), we 
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have the local BH mass density in early-type galaxies: 
P r ly (z = 0) = J M.n c ^ y (M.,z = 0)dM. = / Af.r£ 
= (1.6 ± 0.2) x 10 5 M Q Mpc -3 , 



'((7,0 = 0) 



dcr 

Zm7 



dM. 



(7) 

where a has been equated to a c and da/dAf. is obtained from equation (^). The error estimate represents only the effects of 
the uncertainties in equations (jl|) and (|s]). 



We may also use the mass-luminosity relation i 



1 to estimate p^ arly . 



We note first that the luminosity provided in 



Bernardi 



et al. (2001) is the total luminosity. Early-type galaxies include both elliptical galaxies and SO galaxies. For SO galaxies, the 



bulge and disk both contribute to the total luminosity. Therefore, when applying equation ([]), we need to correct the total 
luminosity of an SO galaxy to its bulge luminosity. Since the BH mass is approximately proportional to the luminosity of the 
hot stellar component (i.e., a.2 = 1.08 in eq. |^ is approximately 1), we may simply correct the total luminosity density of 



early-type galaxies in Bernardi et al. (2001) by the ratio of the luminosity density of hot components to the total luminosity 



density. Fukugita, Hogan & Peebles (1998) estimated the relative luminosity densities fi (in the B band) for elliptical 



galaxies, SO galaxies, spiral and irregular galaxies, and hot components. They found (/e, /so, /spirai+irr) = (0.12,0.23,0.65) 
and /hot = 0.385. They also estimate that in SO galaxies, the bulge luminosity is a fraction /bulge, so//so = 0.64 of the total 
luminosity. Thus, the total luminosity in hot components of early-type galaxies occupies a fraction 

/E + /bulge, so 



_ 0.76 

/e + /so 

of the total luminosity in early-type galaxies. The local BH mass density in early-type galaxies is therefore given by: 

di 



(8) 



P°.t y (* = 0) = 



M.n c ^ ly {M.,z = 0)dM. 



: 0.76 



dM. 



= 2.0 x 10 5 M Q Mpc , (9) 

where the color transformation between g* and B (eq. ^ has been used. 

Massive BHs exist not only in early-type galaxies but also in spiral bulges. Hence, we must augment the estimates (Q) 
and (^) by the contribution from BHs in spiral bulges. As above, we assume that the luminosity fraction is approximately 
equal to the BH mass fraction, and thus the total BH mass density is /hot/(/E + /bulge, so) — 1-44 times the BH mass density 
in early-type galaxies, i.e., from equation (R) 

-3 



p.(z = 0) 



1.44pr rly = (2.3 ± 0.3) x 10 5 M Mpc" 



and from equation 
P.,l{z = 0) 



i.44 P r* 



2.9 x 10 5 M Q Mpc" 



(10) 



(11) 



and (H), we have ignored any 
and m) . Note that equations 



One additional correction to these estimates must be considered. In deriving equations ( 
intrinsic dispersion in the mass-dispersion relation and the mass-luminosity relations (eqs 

(^) and (|^) are fitted in the logM.-log<7 and log M.-log is, hot spaces. Let us assume that the distribution in logAf. at a 
given value of a or Lb, hot is Gaussian, with mean given by the relations (|B|) or (Q), and standard deviation Ai og jv/. that is 
independent of a or Lb, hot- Then our estimates of the mass density must be increased by a factor 



exp 



-(A logA/ . lnlO)' 



l + -(lnl0) 2 A? ogM . + 0(A? ogM . 



(12) 



The intrinsic scatter Ai og m. in the mass-dispersion re lation is small enough that it is difficult to distinguish from measurement 
errors, but is less than 0.27 dex ( Tremaine et al. 2002 ). The corresponding correction factor is therefore between 1 and 1.2, and 
we shall adopt 1.1. The intrinsic dispersion in the mass-luminosity relation is Ai og M. — 0.5 dex (note that only the intrinsic 
dispersion, not measurement error, is included in the definition of Ai og A/. ), so the corresponding correction is a factor of 2.0. 
Thus the density estimates (^) and ([n]) become 



p.(z = 0) = (2.5 ± 0.4) x 10 5 M Mpc" 



and 

P;l(z = 0) = 5. 



x 1O 5 M Mpc" 3 . 



(13) 



(14) 



The error estimate in equation ( |l3| ) represents the combination of the ~10% uncertainty in the correction factor 1.1 due to 
intrinsic scatter, and the error estimate in equation (Jlo|) . The error in equation ( |l4| ) is difficult to estimate but is certainly 
much larger than the error in equation (|13j). Earlier estimates include p*,l(z = 0) = 3.5 x 10 5 Mq Mpc -3 (Salucci et al. 1999), 



P.,l( z = 0) = 3.7 x 10 5 M G 
Salvat|i 200l|). 



Mpc -3 (Merritt fc Ferrarese 2001b), and p., L (z = 0) = (4 ±2) X 10 5 M Mpc~ 3 (Marconi & 



© 0000 RAS, MNRAS 000, 000-000 



Growth of massive black holes 5 



The estimates @ and (Q differ by a factor of more than two. We believe that equation ( |l3| ) is substantially more 
reliable than equation (0), since the correlation between M, and a is tighter and the result does not depend on the uncertain 
bulge-disk decomposition. Nevertheless, the discrepancy between the two equations is worth investigating. The problem cannot 



arise from different estimates of the local galaxy density, since both results employ the same galaxy sample ( Bernardi et al 
20010 nor can it plausibly arise from errors in the BH mass determinations, since the mass-dispersion and mass-luminosity 
relation are based on the same set of mass estimates (see Tremaine et al. 2002 and Kormendy & Gebhardt 2001, respectively). 
To isolate the problem, it is useful to examine the Faber-Jackson relation between dispersion and luminosity. The mean 
relation derived from a large sample of early-type galaxies (Prugniel & Simien 1996) by Forbes & Ponman (1999) is 



log(a/kms _1 ) = 2.25 - 0.102(M B + 20), 



while the mean relation derived from the Bernardi et al. (2001) sample is 
log(o-/kms _1 ) = 2.18 - 0.100(A/ S + 20). 



(15) 



(16) 



Thus the mean absolute magnitude at a given dispersion in the Bernardi et al. (2001) sample is displaced by AiMb = —0.7 
from equation (fl5|). On the other hand, the mean residual in absolute magnitude from the relation ( |l5| ) for the 31 early-type 
galaxies with BH mass determinations in Tremaine et al. (2002) is AjMb = +0.54. A component of the latter difference 
arises because equation (|l^) includes both bulge and disk luminosity in SO galaxies, while the luminosities quoted in Tremaine 



et al. (2002) refer only to the bulges in SO galaxies. Correcting crudely for this using equation (JSp reduces A2MB to +0.24. 
Together, A\Mb and A-^Mb imply that the mean absolute magnitude at a given dispersion is about 1 magnitude brighter 
in the ~ 10 4 SDSS galaxies used to determine the velocity-dispersion and luminosity distribution of nearby galaxies than in 
the ~ 30 galaxies used to determine the mass-luminosity and mass-dispersion relation for BHs. This difference is sufficient to 
explain the difference in BH density estimates in equation s ([l3|) and (fl4|). What is less clear is the origin of this difference. 
Possible explanations include selection effects in either the Bernardi et al. (2001 ) sample or the sample of galaxies with BH 
masses, or systematic differences in the magnitude or dispersion measurements between SDSS and pointed galaxy surveys. 

We close this section with comments on the contribution to the estimated BH density from the faint end of the galaxy 
luminosity function. The SDSS galaxy sample is almost complete at the bright end, but the sample is cut off at the faint end 
(M 9 * — 51og/i > —19.0, cf. Fig. 6 in Bernardi et al. 2001). In our calculation we have extrapolated Bernardi et al.'s distribution 
of magnitude (or velocity dispersion) to arbitrarily faint galaxies. This uncertain extrapolation should not significantly affect 
our estimate of BH mass density, as seen from the following arguments, (i) If we assume that the distribution of absolute 
magnitude obtained in the SDSS is constant for all magnitudes fainter than {M g *} — —19.65 + 51og/i (rather than obeying 
a Gaussian distribution in M s «), the total luminosity density is increased by less than 15%. (ii) The 2dF Galaxy Survey 
extends to considerably fainter absolute magnitudes (Afb, < —14 + 5 log h) [Afbj = M 9t + 0.16 + 0.15(g* — r*) in Norberg et 
al. 2001] than the SDSS sample. The type 1 galaxies in the 2dF Galaxy Survey include all the early-type galaxies and some 
spiral galaxies. The difference between the total luminosity density obtained by assuming that the distribution of absolute 
magnitude of the type 1 galaxies is constant for all magnitudes fainter than Afb T = (M g *) + 0.16 + 0.15(<7* — r*) and that 
obtained by extrapolating the distribution of absolute magnitude of the type 1 galaxies (cf. Madgwick et al. 2001) is less 
than 4%. In addition, in § ^ below (see also Figs. |i| and ^), we will focus on comparing the distribution of local high-mass 
BHs (> 10 8 Mq) with the luminosity function of luminous QSOs, where we need not consider uncertainties of local BH mass 
function at the low-mass ends. 



2.3 BH mass density accreted during bright QSO phases 



Most astronomers believe that QSOs are powered by gas accretion onto BHs, and in this case the luminosity function of 



QSOs as a function of redshift reflects the gas accretion history of local remnant BHs (Boltan 
1992;; ~ 



1982; Chokshi & Turner 



Small & Blandford 1992). We denote the QSO bolometric luminosity produced by a mass accretion rate Af acc as 



Lboi = eAf acc c = eAf.c /(l — e), where e is the mass-to-energy conversion efficiency (up to 31% for thin-disk accretion onto 
a maximally rotating Kerr BH; Thorne 1974), M, = Af acc (l — e) is the growth rate of BH mass, and c is the speed of light. 
The comoving BH mass density accreted during bright QSO phases is given by 



P?, S £(*) = 



tit 
dz' 



dz 



(1 



-*(L,z)dL, 



(17) 



where ^(L,z) is the luminosity function of QSOs, defined so that J ^(L, z)dL gives the comoving number density of QSOs 
at redshift z (L can be either the bolometric luminosity or the luminosity at a given band). In the B band, the bolometric 
correction Cb, defined by Lboi = CbVbL v (B), is about 11.8 (Elvis et al. 1994). Here vbL v (B) is the energy radiated at the 
central frequency of the B band per unit time and logarithmic interval of frequency. 
At z < 3, the QSO luminosity function is often fitted with a double power law: 
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*m(M s ,z) 



10 0.4(fli+l)[M B -M»(^)] + 1Q u 

where ^m(Mb, z)AMb is the comoving number density of QSOs with absolute magnitude in the range [Mb, Mb +dMs] and 



(18) 



*m(M s ,z) = *(L S , z)\dL B /dM B \ = 0.92L B *(L S , z). Boyle et al. (200C) use this functional form to fit the data sets from 
the 2dF QSO Redshift Survey flBoyle et al. 2000| ) and the Large Bright QSO Survey flHewett, Foltz fe Chaffee 1995[ ), which 
contain over 6000 QSOs. For example, in our standard model, (fl m ,f2A) = (0.3,0.7), the luminosity function of QSOs with 
absolute magnitudes —26 < Mb < —23 and redshifts 0.35 < z < 2.3 is described by the following parameters: 



M* B (0) 



*m =2.9 x 10" 5 ft 3 Mpc" 3 mag"\ 
M* B (z) = M* B (0) - 2.5(fci2 + k 2 z 2 ), 
-21.14 + 5 log ft, fci = 1.36, k 2 = -0.27, 
ft = -1.58 and fa = -3.41. 



(19) 
(20) 
(21) 
(22) 

The quadratic dependence of the characteristic magnitude M B (z) on z (eqs. ^ and pi] ) shows an increasing characteristic 
luminosity with increasing redshift at low redshift (z < 2.5) and a decline of the characteristic luminosity at higher redshift, 
which is suggested by observations (e.g. Shaver et al. 1996), although the luminosity function is not yet accurate enough to 
confirm the decline at z > 2.5. The QSO luminosity function over the range 3.6 < z < 5.0 provided in Fan et al. (2001) 
gives a flatter bright-end slope (ft = —2.5) than equation (|li|); however in our calculations below, we will simply extrapolate 
equation @ to high redshift (z > 3) because the detailed QSO luminosity function at z > 3 does not affect our final results 
much (cf. Figure |l] below or the discussion in Chokshi & Turner 1992; if we use the luminosity function in Fan et al. 2001, 
our results change by less than one percent). 

Using equation ( p^ ) we may now obtain the accretion history of the BH mass density arising from optically bright QSO 
phases, p?J°(z). The result is shown in Figure [j], where the mass-to-energy conversion efficiency e is assumed constant. 
About 56% of the present density p^|° (z = 0) is accreted when z < 2 and 90% when z < 3, consistent with our claim that 
P?,acc(z = 0) is insensitive to uncertainties in the QSO luminosity function at high redshift. The accreted BH mass density 
during optically bright QSO phases is found to be 



QSO 



(z = 0) = 2.1 x 10 5 (C* S /11.8)[0.1(1 - ej/e] M Q Mpc" 



(23) 



which is independent of the Hubble constant and remarkably close to the value obtained by Chokshi . 

similar methods, p?,f° (z 
„QSO , 



Turner (1992) using 



0) = 2.2xlO J (CB/16.5)(0.1/e) M Q Mpc" 3 . If the QSOs have a mass-to-energy conversion efficiency 



■ ~ 0.1, p?iuac(z = 0) would be close to the local BH mass density p. (eq. [H|), which implies that the local BH mass density 



comes mainly from accretion during optically bright QSO phases. 

The accreted BH mass density can also be inferred from the X-ray and infrared backgrounds, for example, p; 
9)xl0 5 (0.1/e 



Fabian fc fwasawa (1999[), p. ~ (7.5-16.8) x 10 5 (0 .1/e) M Mpc" 3 in [Elvis, Risaliti fc Zamorani 



;(6- 



Mq Mpc in Haehnelt fc Kauffmann (2001 ) . Multi- waveband (incl uding optical and fa r 



, M Q Mpc" 
(2002fl, p : . R ~ 7.5 x 10 5 (0.1A 
infrared) observations of C han 

|Fabian fc fwasawa (1999 ) conclude that 85% of the emitted energy from AGNs has been absorbed, so that a large fraction 



hard X-ray sources give p. 



HX+OPT+FIR 



9 x 10 5 (0.1/e) M Mpc" 3 in |Barger et al. (200l| ) 



of the local BH mass density must be due to obscured QSOs, which do not contribute to the estimate in equation 



Thus, 



their estimate p^ exceeds p^accO 2 — 0) by a factor 3-4, and then the comparison of p^ with the local BH mass density (eq. |l2 
implies the efficiency e 
Risali 



0.3-0.5, at or beyond the upper limit of plausible accretion processes, e = 0.31 (Thorne 1974) 
Zamorani (2002; ) reached a similar conclusion, e > 0.15, by comparing ~ 
used an early estimate of the local BH mass density 3.5 x 10 5 Mq Mpc 



Elvis 



p;' with the local BH mass density [they 
3 from Salucci et al. (1999), a factor of 1.4 larger than 
our estimate in eq. (|l3|)]. However, there exist many uncertainties in the estimation of the BH mass density from the X-ray 
or infrared background since the population of AGNs contributing to those backgrounds is not yet clearly understood, and is 
not necessarily the same as the optically bright QSOs. 

If the local BH mass density comes mainly from optically bright QSOs, not only the local BH mass density but also the 
local BH mass distribution should be consistent with the accretion history determined from QSOs. The latter comparison is 
the subject of the next section. 



3 SOME CLUES FROM THE BH MASS DISTRIBUTION 

In this section, we employ a continuity equation for the evolution of the BH distribution to relate the local BH mass function 
and the QSO luminosity function. 
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1 2 3 4 5 

z 

Figure 1. The history of the comoving massive BH mass density due to accretion during optically bright QSO phases. The QSO 
luminosity function (eqs. p^|-p2[), which is only well-determined for z<2.5, has been extrapolated to higher redshift. The mass-to-energy 
conversion efficiency e is assumed to be a constant. See details in § £a. 



3.1 Expectations from the continuity equation for BH distributions 

We define S to be any physical variable associated with a BH (e.g. mass, entropy etc.) and ns(S,t) so that ns{S,t)dS is the 
number density of BHs with S — > 5* + dS at time t. The evolution of the distribution ns{S,t) is determined by both BH 
coalescence (caused by galaxy mergers) and other physical processes (e.g. gas accretion) . We can use the following continuity 
equation (cf. Small & Blandford 1992) to describe the evolution of the S-distribution: 

dn s (S,t) d[ns(S,t)(S}} 

1 — 7mcrgcl,0, IJ, 

where 7 mcrgc (S, t) represents the variation of the BH distribution caused by mergers, and (S)(S,t) is the mean rate of change 
of S due to physical processes that conserve the number of BHs. In equation (|2^), we ignore the generation of seed BHs, 
which are usually assumed to have small masses (e.g. < 10 6 M Q ). Since the process of BH mergers is not understood, we shall 
consider two extreme cases: one is the classical case, in which gravitational radiation can be neglected and the total BH mass 
is conserved during a merger; and the other is the adiabatic case, in which gravitational radiation cannot be ignored and the 
total BH area (or entropy) is conserved. We will consider the evolution of the BH mass distribution in the classical case (i.e., 
set S in eq. ^ to be the reducible^ BH mass M.), and the evolution of the area distribution in the adiabatic case (i.e., set S 
in eq. ^i|to be the BH area At). In either case, mergers conserve the total value of S, so we may consider a time interval dt 
in which two BHs parameterized by Si, S2 merge to form a single BH parameterized by Si + S2, i.e., 

7mcr gc (S',t) dt = S(S' -Si- &) - S(S' - Si) - 5(S' - &). (25) 
Thus 

Sf meI ge(S,t) dS = 0. (26) 







If 7mcr gc (S, t) were known, we could integrate equation ( J24| ) to obtain ns(S,t), which could be compared to the local BH 
mass distribution. Unfortunately the BH merger rate is very uncertain. We therefore ask a more restricted question: can we 
find some function of the mass distribution that can only increase during mergers, so that even without knowing the merger 
rate we can establish some inequalities on the current BH mass distribution and QSO luminosity function? For example, we 

1 The dynamically measured BH mass M, (in eqs. jijand ^j) is the reducible mass. For a BH with reducible mass M m and dimen sionless 



spin a (0 ^ \a\ ^ 1), the irreducible mass is given by M ir = M.(l + Vl - a' 2 ) 1 / 2 /%/2 ( Misner, Thorne fc Wheeler 1975| ). Thus, 
M./V2 < M iT ^ M.. 
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would like to choose a function f(S, S') so that 

r meme (S,t) = / /(S,S')7mcrgc(S',t) dS' (27) 



is positive-definite. By inserting equation ( |25[ ) in equation fl27|), we have 

r mcrgc (S,t) dt = f(S,Si + S2)-f{S,Si)-f(S,S 2 ). (28) 

A sufficient co ndition that equation (p8|) is posi tive definite for Si,S2 > is that f(S,S')/S' is a monotonically increasing 
function of S" (Hardy, Littlewood & Polya 1934). We shall choose a specific function f(S, S') below. By integrating equation 
( pi| ) over time t and setting the initial BH S-distribution ns(S,Q) = 0, we may obtain the BH distribution ns(S,to) at the 
present time to as follows: 

^ e [ns (s,tm dt + J\ meigs{s , t)dt (29) 

Now integrate equation ( E?] ) from t = to to, and use equation (|2^) to eliminate 7 mcr gc. Assuming f(S, 0) = 0, we have 

/■*o 

Gioca^S 1 , to) = Guccis, to) + / r mergc (5', t) dt, (30) 

Jo 

where 

/*oo 

Gi c,i(S, to) = f(S,S')n s (S',t ) dS' (31) 



is related to the BH 5-density in local galaxies, and 

G acc (S,t ) = J dS' n s (S',t)(S') df{ ^ S,) dt (32) 
is related to the accretion history of local BHs. Since r mer ge(5', t) is non- negative, we have 

Gi oai i(S,to) ^ Ga, C c(S,to)- (33) 

Inequality (^3|) must hold for every f(S,S') such that f(S,S')/S' is monotonically increasing in S' . We have assumed that 
f(S,0) = 0; it is also useful for f(S,S') to be continuous, since its derivative appears in equation d^). A simple function 
satisfying these conditions, which we shall use for the remainder of this paper, is 

f(S,S') = {°, forS'<S, ( } 

In this case, 

poo 

G loca i(S,t ) = / (S' -S)n s (S',t ) dS', (35) 

poo pt 

Gacc(S,t )= / dS*' / dt. (36) 

Js Jo 

With this definition, equation ( p6[ ) implies that r mcrgc (0,t) = 0, so 

Gi oca i(0,t ) = Gacc(0,to). (37) 
Using this result we can define a normalized version of the fundamental inequality (p3[), 

;c -A ^ — GlocaK-S 1 , to) Gacc(-S', to) _ . . . . 

fflocal(A,t ) = , . ^ , . = gacc(i,t0j- (38) 

(-rlocal(,U, to J Craccl^U, toj 

In the classical case where S = M., equation (^) is identical to Soltan's argument (1982) relating the total local energy 
density in QSOs to the total BH mass density in nearby galaxies, while the inequality fl3q) provides a constraint on the 
normalized BH mass distribution that is independent of Soltan's. Inequalities (^) and ( |38] ) will be our principal tools to 
constrain the accretion history of local BHs. 



3.1.1 Relating the local BH mass function to the QSO luminosity function 

To proceed further, we now assume that the luminosity L during bright QSO phases is an (increasing) function of only the 
central BH mass [dL(M.)/dM, > 0)], which is plausible since usually we have L oc Lboi and, especially for luminous QSOs, 
Lboi(M.) is usually assumed to be near the Eddington luminosity. In our two extreme cases, we have 
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• Classical case: 3 = M, . We have 

G^°(M.,t ) = r dL' P *(L',t) (1 ~ e l LLi dt. (39) 

Jh{M.) Jo eC 

which represents the time-integrated mass density accreted onto BHs with masses larger than M m during optically bright QSO 
phases. Equation (^) is independent of the duty cycle of bright QSO phases, so long as the luminosity during these phases 
is determined by M m . If the local BH mass density arises entirely by mergers and by accretion during optically bright QSO 
phases, then by combining inequality (|3^) with equation ( pS9| ) we have 

Gg? c °(M.,i ) = Gacc(M.,to), (40) 

mocal(M "' 0) = G local (0,i ) > GT°(0,t ) = 9 - (M " t0) - (41) 

• Adiabatic case: S = A, — WirG 2 M?/c 4 , where G is the gravitational constant and Mi r is the irreducible mass of a BH 
with reducible mass M. (see footnote 1). We assume that the dimensionless spin of BHs in QSOs is in an equilibrium state 
(i.e., a — 0), which we denote by (Jqso- For example, in the thin-disk accretion model for QSOs the dimensionless spin a 
increases from to reach the equilibrium value ciqso = 0.998 after the BH accretes ~ 1.5 times its initial mass (Thorne 1974). 
Thus, in the equilibrium spin state, we have the area (or entropy) variation rate of a BH as follows: 



A. = 32nG 2 M il M il /c 4 = 16tt G 2 ' M.M. (l + Jl- a% so ) /c (a = 0). (42) 

If the local BH entropy density arises entirely by mergers and by accretion during optically bright QSO phases, then equation 
( p| ) can be written as 

G m (A.,to) = dL' 9(L',t)A.(M'.,a Q so) dt, (43) 

J C(M ir ) JO 

where the function £(Mi T ) is defined by £(Mi r ) = L(M.). It proves useful to use the irreducible mass rather than the area as 
the independent variable, so we define 

gtc' QS °(M iT ,t ) = G acc (A.,t ), where A. = 16ttG 2 M 2 /c 4 . (44) 
In the adiabatic case, equation ( pB[ ) may be written as 

l> oo 

Giocai^Wo) = / [A.{M' m ,a)-A.\n M .{M' m ,t )dM' m (45) 

J M.(A.,a) 

where for the moment we assume that all the local BHs have the same spin a. Unfortunately, the BH spin parameter a in 
dead quasars is unknown (in particular, there is no reason to suppose that it equals the equilibrium spin parameter oqso for 
accreting BHs). However, for a given value of the reducible mass M, , both the irreducible mass Mi r and the area A. decrease 
with increasing a. Thus we have the inequality 

Gi oca i(A.,to) ^ / [A.(M'„a = Q)~ A.]n M .{M'.,t )dM'.. (46) 

Jm.{A.,o) 

Once again introducing the irreducible mass rather than the area as the independent variable, we define 

Qi'cAMivM) = / {A.(M'.,a = 0)-A.(M lI )]n M .(M'.,t )dM'.. (47) 

J Mi! 

The integral on the right side of equation ( [fiil) differs from the integral on the right side of equation ( [f7j) only in the lower 
limit of the integrand. Since M[ T ^ M, (see footnote 1), the second integral has a larger integration range and therefore is 
larger. Thus 

G local (^.,t ) SS 01ocal( M ir,to). ( 48 ) 

The inequality ( fii| ) is an equality if and only if BHs in local galaxies are all Schwarzschild BHs (a = 0). By using inequalities 
(p|) and (E|) and equations (pi) and (pi), we have finally: 



2 We assume that any QSO is due to accretion onto a single massive BH, rather than accretion onto two or more BHs in the galactic 
nucleus, although binary BHs are likely to be present in some galactic centers (Begelman, Blandford & Rees 1980; Yu 2002). If there 
exists more than one BH in a QSO and each BH is accreting with the Eddington luminosity, the mass M m obtained from the total QSO 
luminosity will represent the total mass of the BHs, and our arguments in the classical case are unaffected. However, in the adiabatic 
case, the entropy or area increase associated with a given mass accretion rate depends on the number of BHs. 
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GtAMirM > a a A c -' QS °(M ir , to). (49) 



Note that inequality (J49j) still holds even if we discard the assumption made after equation ( |45[ ) that all local BHs have the 
same spin parameter. 

To estimate Gi oca i and £? lo * al , we proceed as follows. High- luminosity galaxies are mainly early- type galaxies (see, for 
example, Figure 4.14 in Binney & Merrifield 1998 or Figure 7 in Madgwick et al. 2001); moreover, the mass of the central 
BH is approximately proportional to the luminosity of the hot stellar component (cf. eq. []), which is less than ~30% of the 
total luminosity for spiral galaxies (Fukugita, Hogan & Peebles 1998); hence high-mass BHs (M. > 1O 8 M0, corresponding 
to MB,hot < —20) are mostly in early-type galaxies. The local BH mass distribution in early-type galaxies can be obtained 
from their velocity dispersions and the mass-dispersion correlation. By analogy with equations Q35|) and (|45j), we define 



C(W.,!o)= (M:-M.)n^(Mi,t )dMi (50) 

J M, 

and 

Gtr ly (M iT ,t ) = / [A.(Mi,a = 0) - A.(M ir )K^ ly (i\Oo) AM 1 .. (51) 

J M. 

For high-mass BHs (M. > 10 8 Mq), we have GZ^(M.,t ) ~ Gi oca i(M„ t ) and 0£j£ arly (M.,f Q ) ~ G; A o ' c JM.,t ). In 
addition, as seen from inequalities (§|), @ and @, we have 

G^(M.,to) > GT°(M.,to), (52) 

early/ , \ ^local 

ffw (M.,t ) = Giocai(0i , o) £ G QSO (0ito) = ( M -*°) ~ 10 Me) ' (53) 
in the classical case (with approximate equality if and only if mergers are neglected) and 

Otr ly (Mir, t ) > g£v QS °(M iT , t ) (M ir > 10 8 M ) (54) 

in the adiabatic case (with approximate equality if and only if mergers are neglected and BHs in local galaxies are all 
Schwarzschild BHs). The approximate nature of the inequalities arises solely from our restriction to early- type galaxies. 



3.2 Comparison with observations 

We now compare inequalities (^)-(^) with the observations. Throughout this section, all the variables related to the dis- 
tribution of BH mass and entropy density, G^(M.,t ), g^cl ( M *^o) and G^'cS^ ( M ^, to), are obtained by using the 
mass-dispersion relation (eq. []) and ignoring the small intrinsic scatter in this relation. 



3.2.1 The classical case 

In the classical case, the total BH mass is conserved during mergers, which is implicitly assumed in Soltan's argument and 
many subsequent discussions of the growth of BHs. In § [| we have obtained that the total BH mass density in local galaxies is 
consistent with that accreted during optically bright QSO phases, if the mass-to-energy efficiency e ~ 0.1, i.e., Gi OC ai(0, to) ~ 
Gj3?°(0, to) (cf. eqs. ^ and ^) . In contrast, inequality ( p^ involves normalized quantities g(M,,to) = G(M m ,to)/G(Q, to) and 
therefore provides a constraint independent of Soltan's argument. [For example, note that the normalized quantity g acc (M,,to) 
in equation (53) is independent of the efficiency e if e is a constant.] 

Figure |2 (a) shows G^^(M. , to) derived from the SDSS sample of early-type galaxies (dotted and dashed lines, eq. |5c| ) 
and G^ C °(M. , to) derived from optically bright QSOs (thick solid line, eq. |3{]). The QSOs are assumed to have the Eddington 
luminosity and the efficiency e is given by e/(l — e) = 0.1. Figure ^|(b) shows the normalized quantity g^^iM,, to) derived 
from the SDSS sample of early-type galaxies (dotted and dashed lines, eq. ^) and g2^°(M,, to) derived from optically bright 
QSOs (thick solid line, eq. ^) using the same assumptions. The dotted line in panel (a) [or panel (b)] is used for G^*^(M,,to) 
[or ff^cal (-^j to)] in the region M. < 10 s Mq where spiral bulges can make a substantial contribution, so that G^J (or g^*^) 
is only a lower bound to Giocai (or giocai). For G^'^(Af., to) and g\^J^{Mm,to), we show separately the results obtained by 



using the mass-dispersion relations i n [Tremaine et al. (2002 ) (thick dotted and dashed line) and Merritt & Ferrarese (2001a) 
(thin dotted and dashed line). Figure g shows that for high BH masses (M. > 10 s Mq), G^ c °(M.,t ) [or g^ s c ° (M.,t )] (thick 
solid line) is larger than G^!^(A/. , tp) [or flf^j to)] (dashed lines) by a factor of more than two, which is inconsistent 
with inequalities (^) (or ^3|) and (^) . 

This inconsistency is robust to a number of changes in our assumptions. For example: 
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Figure 2. Panel (a): Variables related to the BH mass density as a function of BH mass M m . The thick solid line shows G? C c°(iW., t ) 
contributed by QSOs (eq. [ji) ), where QSOs are assumed to have the Eddington luminosity and the mass-to-energy efficiency e is given 
by e/(l - e) = O.f . The dashed line (M, > fO 8 M ) and the dotted line (M. < 10 8 M Q ) show G™£* (M.,t ) obtained from the SDSS 
sample of nearby early-type galaxies (eq. pc| ). In the dotted region G^^l un d eres ti ma tes G\ oca \ because of the contribution from BHs 



in spiral bulges. The thick dotted and dashed lines show the result obtaine d by using the mass-dispersi on relation in Tremaine et al 
(2002j)J and the thin lines show the result obtained by using the relation in Merritt 

to the normalized BH mass density as a function of BH mass Mm, i.e., g^x (Mm, to) and 9,°""']' to) (eq. |53| 



(2001a| ). Panel (b): Variables related 
QSOs arc assumed to 



>*°) and ~9l^f( M «'* ) ( ec l _ 

have the Eddington luminosity and the curves are independent of the mass-to-energy efficiency e, so long as it is a constant. The line 
t ypes have the same meanings as those in Panel (a). At the low-mass end, g^^(M,,to) approaches (/e + /bulge, so)//hot — 0.69 (see 
§ 2.2). For high BH masses (M. > 1O 8 M ), inequalities (||) and (||) imply that the thick solid line should be equal to the dashed 
one in panels (a) and (b) if mergers are unimportant, and lower if mergers are important, which is inconsistent with the observations. 
The thin solid line represents an estimate of G^ C c°(M. , to) in panel (a) [or g2c° (Mm , to) m panel (b)] obtained by cutting off the QSO 
luminosity function at Lbol J/ f.hh(10 9 Mm), and th is est imate is still significantly higher than the dashed lines and inconsistent with 



- — - ~..-~. J -~— — DUl - ~DUU y 

inequality (B2J) (or inequality |53[) . See discussion 



3.2.1 
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• We have assumed that QSOs radiate at the Eddington limit; if they have sub-Eddington luminosity (e.g. one tenth of 
the Eddington luminosity in Ciotti, Haiman & Ostriker 2001), the solid lines in Figures ^(a) and (b) would shift to higher 
mass and the discrepancy would become worse. 

• The QSO luminosity function is uncertain at the bright end, and the mass-dispersion relation is uncertain at the high- 
mass end. However, even if we cut off the QSO luminosity function at Lboi > L E dd(10 9 M ) ~ 10 47 erg s -1 (L Edd : Eddington 
luminosity), the new G^. (M.,t Q ) [or g^ c s c °(M.,t )] (thin solid line in Figure ||) is still significantly higher than (M., t ) 
K <C£(M.,t )] for 10 s M < M. < 10 9 M . 

• If there exist mechanisms other than accretion during bright QSO phases [e.g. accretion via advection dominated accretion 
flow (ADAF), see discussion in Haehnelt, Natarajan & Rees (1998); or accretion of non-baryonic dark matter] that contribute 
to the local BH mass density, then G^ C |°(M. ,to) is only a lower limit to G aC c(M. ,t), and the inconsistency with inequality 
( p3[ ) for high-mass BHs would become worse. 

• Study of the X-ray background suggests that a majority of the total BH mass density is contributed by obscured accretion 



(Fabian & Iwasawa 1998). If the number ratio of obscured QSOs to optically bright QSOs is independent of BH mass (e.g., 
in the standard unification model for AGNs, obscuration is a purely geometrical effect and obscured fraction has no relation 
to central BH mass), gr a cc(M. ,to) would be unchanged, and the inconsistency in Figure ^(b) would remain. 

One possible resolution to the inconsistency shown in Figure ^ is that the mass-to-energy conversion efficiency e depends 
on QSO luminosity. In particular, the average efficiency in luminous QSOs (e.g. Lboi > 10 46 ergs -1 , which corresponds to the 
Eddington luminosities of BHs with M, > 10 s Mq) must be > 0.2 to make the BH mass density accreted during optically 
bright QSO phases as large as that in local galaxies (cf. Fig. ^a). Note that the required average efficiency is close to the 
maximum efficiency (~ 0.3 in Thorne 1974) allowed in thin-disk accretion models. This result suggests that other accretion 
mechanisms arc less signifi cant than accretion during optically bright QSO phases, and we may expect that obscured accretion 



( Fabian fc Iwasawa 1999 ) suggested by the X-ray background is not important for the growth of high-mass BHs, which is 
consistent with the scarcity of Type II QSOs. 

If the growth of low-mass BHs also occurs mainly during optically bright QSO phases, then less luminous QSOs (-Lboi ^ 
10 46 ergs -1 ) must accrete with an efficiency less than ~ 0.1 to maintain the consistency between the local BH mass density 
and the BH mass density accreted during b right QSO phases (cf. eq. |3^ ). If on the other hand, obscured accretion contributes 
significantly to the total BH mass density ( Fabian fc Iwasawa 199s| ), less luminous QSOs may have efficiency ~ 0.1 or even 



greater. In this case, the obscured fraction must be larger for low-mass BHs than for high-mass BHs, which might suggest 



that obscured accretion occurs mainly for small BHs in young growing galaxies as described in the model by Fabian (1999). 
Other possible resolutions to the inconsistency shown in Figure ^ include: 

• Not all massi ve BHs ma y reside in galactic centers; for at least three reasons: (i) BHs may be left in the halo after galaxy 
mergers; however, Yu (2002; ) shows that the orbital decay time from dynamical friction for BHs in the galactic halo is generally 



much less than a Hubble time, so long as the orbiting BHs remain part of the tidally stripped remnant of their original host 
galaxy, (ii) Incomplete merging of binary BHs formed in mergers is expected to result in binary BHs with semimajor axes in 



the range 10 -3 -10 pc depending on the velocity dispersion and shape of the host galaxy, and the masses of the BHs ( Yu 2002); 
however, in most cases the separation is so small that the binary should look like a single merged BH at current telescope 
resolutions, (iii) BHs can be ejected from galaxy centers through either three-body interactions with other massive BHs (e.g. 
Valtonen 1996) or gravitational radiation reaction during BH coalescence (e.g. Rees 2001). 

• The bolometric correction that we have used, Cb = 11.8, may be too large; however, the uncertainty in this estimate of 



the bolometric correction (Elvis et al. 1994) is only ±4.3, and the smallest bolometric correction in the Elvis et al. sample is 



5.5, so it is unlikely that errors in the bolometric correction can lower the required efficiency by more than a factor of two. 



Luminous QSOs may accrete with super- Eddington luminosities (Begelman 2001, 2002) 



3.2.2 The adiabatic case 

BH mass may not be conserved during BH mergers, since mass can be radiated as gravitational waves. According to the 
second law of BH thermodynamics, BH mergers never decrease total BH entropy. Here, we study an extreme case in which 
the BH entropy is conserved during BH mergers. 

Figure^] shows the variables related to the BH entropy density as a function of BH irreducible mass that we introduced 
in § |3.l| : G lo l£ arly (Mi r ,to) obtained from the SDSS sample of early-type galaxies (dotted and dashed lines, see eq. ^) and 
G^ c ' c ' c>so (Mi I ,to) obtained from optically bright QSOs (thick solid lines, see eq. ^). As usual QSOs are assumed to have the 
Eddington luminosity and e is given by e/(l — e) = 0.1 (eqs. |li|-p2|). The spin of BHs in QSOs is assumed to be in an equilibrium 
state, a — 0, a = aqso- The upper and lower thick solid lines represent the results obtained by assuming that all the BHs in 
QSOs are Schwarzschild (aQso = 0) and maximally rotating Kerr BHs (aQso = 1), respectively. The dotted lines are used for 
5 1 ' 4 *^j arly (A/i r , to) in the region Mix < 10 8 Mq where spiral bulges can make a substantial contribution, so that S 1 ^* a ° arly is only 
a lower bound to 5^* al - For Qi l^ ly (Mix, to), we show separately the results obtained by using the mass-dispersion relations 



© 0000 RAS, MNRAS 000, 000-000 



Growth of massive black holes 13 



Tremaine et al. (2002) (thick dotted and dashed line) and in Merritt fc Ferrarese (2001a| ) (thin dotted and dashed line). 
The functions G^^ 1 ^ (M- lr , to) obtained from these two relations differ by a factor of two, much more than the analogous 



curves in Figure |2j. The difference is mainly caused by BHs with M. > 10 9 Mq because the BH area (or entropy) A. oc M? so 
the results are quite sensitive to the mass-dispersion relation at the high-mass end. Figure ^ shows that the thick solid lines 
[i.e., G&'c'^ {M\ T , to)] are much higher than the dashed lines [i.e., ^*^i arly (A^ir, to)], no matter whether we choose oqso = 
or aQso = 1 and for both mass-dispersion relations, which is inconsistent with inequality (Q). This inconsistency may be 
partly caused by uncertainties in the QSO luminosity function at the bright end. To investigate the effects of uncertainties in 
the QSO luminosity function, we cut off the luminosity function at Lboi > LEdd(10 9 Mq) (thin solid lines in Figure [|). For 
consistency, we must then use the Tremaine et al. mass-dispersion relation since it gives £?i o £ii arly (10 9 M0,to) — 0. We see 
from Figure ^ that the new <5^c' < ^ SO (AL. r , to) for aQso = (upper thin solid line) is still significantly higher than the thick 
dashed line, which is inconsistent with inequality (|54|). Considering that the maximum efficiency for thin-disk accretion onto 
Schwarzschild BHs is less than 0.1, a more realistic S^ ,QSO (M ir , t ) for iqso = (i.e., e ~ 0.06) should be even higher than 
that shown in Figure |§| and the discrepancy between the upper thin solid line and the thick dashed line would become worse. 
However, the new 5^J' CiSO (Mi r , to) for ceqso = 1 (lower thin solid line) is close to or only slightly higher than the thick dashed 
line, which is not seriously inconsistent with inequality (pi|), and any inconsistency can be reduced further because rapidly 
rotating BHs can have efficiencies substantially larger than 0.1. Thus, Figure ^ suggests that: it is impossible that all BHs 
in luminous QSOs (Lboi > 10 46 ergs" 1 ) are Schwarzschild BHs; BHs in most luminous QSOs are rapidly rotating Kerr BHs; 
and luminous QSOs are accreting at an efficiency > 0.1 (cf. Elvis, Risaliti & Zamorani 2002). 

Possible solutions to the inc onsiste ncy (mainly for the comparison of local BHs with QSOs having <iqso = 0) shown in 
Figure [] have been discussed in § 3.2.1 An additional possibility is that there exists more than one BH accreting materials in 
a QSO. In this case, the entropy or area increase associated with a given mass accretion rate depends on the number of BHs 
(see footnote 2) and equation ( |42| ) is an overestimate of the area increase. 

In the adiabatic case, we do not discuss the roles of obscured accretion (or other accretion mechanisms) on the BH entropy 
density since the total BH entropy obtained from the X-ray background is not clear and also the conclusion may be affected 
by the uncertainty of the mass-dispersion relations at the high BH-mass end or the QSO luminosity function at the bright 
end. 



3.3 Mass-to-energy conversion efficiency and mean lifetime of QSOs 



In S 3.2 



we have argued that the growth of BHs with M. > 10 Mq mainly occurs during optically bright QSO phases and 
the mass-to-energy conversion efficiency of luminous QSOs should be larger than 0.1. In this subsection, we will ignore BH 
mergers [the term Ymerge(M. , t) in eq. ^] and estimate in more detail the mass-dependent efficiency of QSOs that is required 
to match the local BH mass function. We will also give an estimate of the mean lifetime of QSOs. 



3.3.1 Mass-to-energy conversion efficiency 



We assume that the mass-to-energy conversion efficiency is a function of BH mass e(M.), or a function of QSO bolometric 
luminosity e(Lboi) since Lboi is assumed to be a function of only BH mass M. . By replacing the variable S in equation ( |24[ ) 
by M' m and ignoring the term 7 mGr gc, we have 



dn M .(M'.,t) d 



dt 



+ 



DM'. 



9(L,t) 



AL (1 - e)L hol (M'.) 
AM'. ec 2 



0. 



By integrating equation (^5|) over M' m from M. to oo and over t from to to, we find 
e(Lboi) = 7/(Lboi)/[l + r;(Lboi)], 



(—) 
\AM.J 



J t0 *(L,t) At 



(55) 



(56) 



(57) 



, and assuming 

that QSOs have a fixed fraction of the Eddington luminosity (we examine two cases, 0.1 and 1), we obtain the efficiency e(-Lboi) 
in Figure ^. The solid lines in the interval 10 46 ergs _1 < Lboi < 10 47 ergs _1 [corresponding to the BH mass in the range 
10 s (L Bd d /Lboi )M < M. < lO 9 (L E dd/L b oi)M ] denote the region where the observational constraints are fairly strong. The 
dotted lines at the luminous and faint ends are obtained by extrapolating the BH distribution in early-type galaxies or the 
QSO luminosity function. We also show separately the results obtained by using the mass-dispersion relations in Tremaine et 



where 



J?(L b oi) 



n{M'.,to)AM'. 



r oo 
JM.(L bol ) 

By inserting the QSO luminosity function and the BH mass function in local early-type galaxies in equation (|56| 



al. (2002) (thick dotted and solid lines) and in Merritt & Ferrarese (2001a) (thin dotted and solid lines). The lower thick/thin 
solid and dotted lines represent the results by assuming Lboi = L-^dd', and the upper ones represent the results by assuming 
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Lboi = O.ll/Edd- The fact that the efficiency is required to be higher than 0.3 (the maximum efficiency for thin-disk accretion 
onto BHs) or even approach 1 when Lboi = O.lLEdd suggests that the assumption of sub-Eddington luminosity for luminous 
QSOs is unrealistic, a conclusion derived by different arguments in § 3.2.1. Figure ^] shows that even for Lboi = LEdd, when Lboi 



approaches 10 ergs - the mass-to-energy efficiency may be slightly larger than 0.3 (the maximum efficiency for thin-disk 
accretion onto BHs). As described in § 3.2.1, this unrealistic e may be caused by (i) ignoring BH mergers, (ii) ignoring ejection 
of BHs from galactic nuclei or BHs left in the halo after a galaxy merger, (iii) uncertainties in the mass-dispersion relation 
at the high-mass end or the QSO luminosity function at the bright end, or (iv) ignoring the possibility of super-Eddington 
luminosities. 



3.3.2 Mean lifetime 

If the duty cycle of QSOs with BH mass M. is 5(M.,t) at time t, we may define the mean lifetime of QSOs as follows: 

t(> M.) = ^ . (58) 

With this definition, we have 

1 J~n M .(Mi,t )dMi - f^n™*(M^,t )dMl ' ^ 

the approximate equality holds for M, > 1O 8 M0, where early- type galaxies dominate the number density. The lifetime 
t(> M.) in equation (^9|) is independent of any assumed value of the efficiency e. By inserting the QSO luminosity function 
and the BH mass function in local early-type galaxies in equation (|5!]), and assuming that QSOs have the Eddington luminosity, 
we obtain the lifetime r(> Al,) as a function of M, in Figure [| As seen from Figure ^, the mean lifetime obtained by using 



the mass-dispersion relation in Tremaine et al. (2002| ) (thick solid line) is r(> M.) ~ (3-13) x 10 7 yr, i.e., the average duty 



cycle of QSOs is about r/tnubbic — 3-13 xlO -3 , where tHubbic is the Hubble time ~ 10 10 yr. This range is comparable to the 
Salpeter time ~ 4.5 x 10 7 e/[0.1(l — e)] yr, the time for a BH accreting with the Eddington luminosity to e-fold in mass (for 
a BH with spin in an equilibrium state, the time for a BH radiating at the Eddington luminosity to e-fold in entropy is half 
of the Salpeter time). Hence, accretion during bright QSO phases can significantly increase the mass of BHs, which is once 
again consistent with the argument that growth of high-mass BHs comes mainly from accretion during optically bright QSO 
phases. The lifetimes of QSOs derived here are also consistent with some theoretical models for the luminosity function of 



QSOs (e.g. 


iaehnelt, Natarajan & Rees 1998 


& Hui 2001 


j 


Martini & Weinberg 2001 


)■ 



Haiman & Loeb 1995 ) and with the spatial clustering of QSOs (e.g. 



Haiman 



4 CONCLUSIONS 

In this paper, we have studied the observational constraints on the growth of massive BHs in galactic nuclei. We use the 
velocity dispersions of early-type galaxies obtained by the SDSS and the relation between BH mass and velocity dispersion 
to estimate the local BH mass density to be p.(z = 0) ~ (2.5 ± 0.4) x 10 5 /ig. 65 M Mpc" 3 (eq. [t|). We also use the QSO 
luminosity function from the 2dF QSO Redshift Survey to estimate the BH mass density accreted during optically bright 
QSO phases to be (z = 0) = 2.1 x 10 5 (C S /11.8)[0.1(1 - e)/e] M Mpc" 3 (eq. ||). These two results are consistent if the 
QSO mass-to-energy conversion efficiency e is ~ 0.1. 

By studying a continuity equation for the BH distribution and including the effect of BH mergers in two extreme cases 
(the classical case and the adiabatic case), we predict relations between the local BH mass function and the QSO luminosity 
function. In the classical case, the predicted relation is not consistent with the observations at high BH masses (M. > 10 s M ) 
unless luminous QSOs (Lboi > 10 46 erg s" 1 ) have an efficiency higher than 0.1 (e.g. e ~0.2, which is possible for thin-disk 
accretion onto a Kerr BH; see Elvis, Risaliti & Zamorani 2002 for a similar conclusion). Possible ways to evade this conclusion 



include super-Eddington luminosities (Begelman 2001, 2002) or ejection of BHs from galactic centers, but if these are not 



important then we will come to the following conclusions. Accretion other than what is traced by optically bright QSOs is not 
important for the growth of high- mass (> 10 8 M©) BHs (e.g. accretion by obscured QSOs, accretion via advection dominated 
accretion flow, or accretion of non-baryonic dark matter). If the growth of low- mass BHs also occurs mainly during optically 
bright QSO phases, less luminous QSOs (Lboi < 10 46 erg s _1 ) must accrete with a low efficiency < 0.1; alternatively, other 
mechanisms such as obscured accretion may contribute to the mass density of low-mass BHs, and their efficiency could be 
similar to high-mass BHs, e ~0.2. The comparison of the observations with the predicted relation in the adiabatic case also 
suggests that BHs in most luminous QSOs should be rapidly rotating Kerr BHs with an efficiency > 0.1, or there should 
exist more than one BH in a QSO. We may expect that these constraints on BH demography are useful for estimating the 
gravitational wave signal from BH mergers that may be detected by Laser Interferometer Space Antenna (LISA). 
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We estimate the mean lifetime of luminous QSOs with BH mass in the range 10 s M < M. < 10 9 M Q , which is (3- 
13)xl0 7 yr and comparable to the Salpeter time if e ~0. 1-0.3, which is once again consistent with the argument that the 
growth of high- mass BHs {M. > 1O 8 M0) comes mainly from accretion during optically bright QSO phases. 

We are grateful to Amy Barger, Mitchell Begelman, Xiaohui Fan, Jeremy Goodman, Martin Haehnelt, Zoltan Haiman, 
Youjun Lu, David Schlegel, David Spergel and Michael Strauss for helpful discussions. This research was supported in part 
by NSF grant AST-9900316 and by NASA grant HST-GO-09107.09-A. 
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Figure 3. Variables related to the BH entropy density as a function of BH irreducible mass M\ v , i.e. 



,to) (eq. 0, 



thick 



solid lines) and Q^ > '^ lly (Mi I ,to) (eq. |5l| dotted and dashed lines). QSOs are assumed to have the Eddington luminosity and an efficiency 
e given by e/(l — e) = 0.1. The dotted/dashed lines and the thin/thick lines have similar meanings to those in Figure H. The spin of BHs 
in QSOs is assumed to be in an equilibrium state ciqso- The upper thin/thick solid lines are obtained by assuming that all the BHs in 
QSOs are Schwarzschild BHs (aQgo = 0) and the lower thin/thick solid lines assume that the BHs are maximally rotating (ciqso = !)• 
For high BH masses (M^ > 10 8 Mq), inequality ( pi} ) implies that the thick solid line should not be higher than the dashed one. After 
cutting off the QSO luminosity function at Lbol ?J ^Edd(10 9 Mq), the upper thin solid line for ciqso = is stil l significantly higher than 
the thick dashed line, which implies that the BHs in QSOs must be rapidly rotating. See discussion in § 3.2.2 
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Figure 4. The (bolometric) luminosity-dependent efficiency of QSOs e(Lbol) if mergers are unimportant (eq. ^). The solid lines in the 
range 10 46 ergs -1 < Lbol ^ 10 47 ergs -1 [corresponding to the BH mass in the range 10 8 (-LEdd/^bol) Mq < M, < 10 9 (LE ( jd/i'bol) Mq] 
show the region in which the observational constraints are fairly secure, and the dotted lines outside this region show the results obtained 
by extrapolating the BH mass function in local galaxies and th e QSO luminosity function. The thick dotted and solid lines are obtained 
bv using the mass-dispersion relation in Tremaine et al. (2002), and the thin dotted and solid lines are obtained by using the relation in 
tvlerritt & Ferrarese (2001a). The lower thick/thin lines represent the results by assuming Lbol = Lgdd, and the upper ones represent 
the results by assuming Lbol = O.lLEdd- If Lbol = O.IL waa. t he efficiency is required to be higher than 0.3 (the maximum efficiency for 
thin-disk accretion onto BHs) or even approach 1. See 



3.3.1 for discussion. 
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Figure 5. The mean lifetime of QSQs with BH mass > M, as defined by equation (|59|), The thick/thin and solid/dotted lines have 
similar meanings to those in Figure W. The horizontal dashed lines repr esent the Salpeter time (the time for a BH radiating at the 
Eddington luminosity to e-fold in mass) ~ 4.5 X 10 7 e/[0.1(l — e)] yr. See § 3.3.2 for discussion. 
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